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Abstract 

We extend an earlier model for radiatively generated fermion masses based 
on the Pati-Salam group to include CP violation. Spontaneous CP violation in 
the early universe gives rise to a complex mass matrix for heavy sterile neutrinos. 
The out-of-equilibrium decay of these neutrinos generates a B — L asymmetry. 
The sterile neutrinos also act as a mass seed in generating one-loop (complex) 
mass matrices for the quarks. Thus, the two low energy manifestations of CP 
violation - the CKM phase and the baryon number asymmetry - can both be 
traced in a calculable way to a common source. 
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1 Introduction. 



All of the information we have about CP violation can be summarized by two 
seemingly unrelated numbers: e = 2.26 ± 0.02 x 10~ 3 and rjs = ns/n^ = (2.8 — 
4.0) x 10 -10 |]|. These two numbers are a measure of two seemingly unrelated 
phenomena: the strange decays of the K mesons and the baryon asymmetry of 
the universe. It is generally agreed upon that the best candidate for the source 
of e is an unremovable phase in the quark mixing (CKM) matrix ||; however, 
there is no such agreement on a single most likely source of t]b- 

Although there are a number of viable models of baryogenesis most of them 
do not say anything about the CKM phase, a notable exception to this being 
some of the papers that propose that the baryon asymmetry is generated at the 
electroweak phase transition with the only source of CP violation being in the 
quark mass matrix ||. Although electroweak baryogenesis remains an exciting 
possibility 0] the general consensus is that additional sources of CP violation are 
needed to make it work ||. 

An interesting mechanism of generating the baryon asymmetry is to first 
generate a lepton asymmetry by the out-of-equilibrium decay of a heavy sterile 
neutrino and then use the B + L anomaly in the Standard Model to cycle this into 
a baryon asymmetry || [7] . This scenario fits naturally into a model for radiatively 
generated fermion masses based on the Pati-Salam gauge group SU (4) x SU (2) l x 
SU(2)r II that we have presented ||. The motivation for this model was firstly 
to treat the quarks and leptons on a symmetric footing, and then to show that 
one can generate realistic masses and mixings for the Standard Model fermions 
without a large hierarchy of Yukawa couplings. 

In this model the Standard Model fermions (including right-handed partners 
for the neutrinos) are supplemented by three generations of heavy sterile neu- 
trinos. As a result of a particularly simple choice of scalar representations it is 
impossible for the Standard Model fermions to get masses at tree-level, and the 
sterile neutrinos act as mass seeds in generating finite one-loop masses for them. 
If we allow the sterile neutrinos to have a complex mass matrix (as the result of a 
spontaneous violation of CP in the early universe), then not only does their out- 
of-equilibrium decay generate a baryon number asymmetry, but their exchange 
in the one-loop diagrams also generates a complex mass matrix for the quarks. 
Thus one is able to derive both e and t]b from the same source: a complex phase 
in the sterile neutrino mass matrix that gets communicated to the quark sector 
at the scale where lepton number breaks down. It is the purpose of this paper to 
explore this possibility. 

In a sense this model is complementary to the models of electroweak baryo- 
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genesis where one starts with the CKM phase and tries to use it to generate 
a fermion number asymmetry. Here we start with a fermion number violating 
operator and try to use it to generate a phase in the CKM matrix. 

In Sec. 2 we briefly review the model, its representation content and pattern 
of symmetry breakdown. In Sec. 3 we calculate the baryon asymmetry generated 
in this model in two alternative scenarios: one in which there was a period of 
inflation in the early universe, and one in which there wasn't. In Sec. 4 we discuss 
the fermion masses and CKM matrix, and illustrate with a numerical example. 
We conclude in Sec. 5. Some details of the scalar potential and explicit formulas 
for the fermion masses are presented in Appendices A and B. 



2 The Model 

The gauge group is S77(4) x SU(2) L x SU{2) R with gauge couplings g$, 9l and 
gji. The Standard Model fermions transform in the usual representations: 

n ~ (4,2,1)' s ( l l 2 e v I (i) 

(4,1,2)., (I 2 I »)[ (2) 

where i — 1,2, 3 is a generation index, and we have included a right handed 
neutrino N. We add to this three generations of (right-handed) sterile neutrinos 

**~ (1,1, 1)*- (3) 

We choose to make the matter spectrum supersymmetric, so the scalars in the 
model transform as 

L< ~ (4,2,1)' = ( L T ul L T u2 f" 3 L T V \ , (4) 



Rui Ru2 Ru3 Rn 



i 



K ~ (4, 1,2)»= ^ ^ T , (5) 

V n dl -Kd2 J^d3 n e ) 

and 

^ ~ (1,1,1)'- (6) 



1 This differs from the scalar spectrum in Ref. where we had only two generations of L 
and R, and no a. For other related models see Refs. H, 0. 
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We will impose a discrete Z% symmetry on the gauge singlets (broken by the 
interactions of the Standard Model particles) under which s J — > e~ l ^^^s^ and 
This permits us to make the Lagrangian CP invariant, with the 



vacuum expectation values of the a 3 breaking CP spontaneously as in Ref. |T2 
The Yukawa interactions will then be 

Cy = -y t (s c YsW - K)#is^° - {KnW^R* + h.c, (7) 

with all of the coupling constants real. We discuss the details of the scalar 
potential in Appendix A, noting only that we can choose parameters such that 
it is minimized when 



v ° j*i I p \ Vr x I r \ Vl ; 
-j=e J , {K N ) j = —dy, {L u )j = -j=c 



with |uq| > \vr\ ^> \vl\- None of the Standard Model fermions get masses at tree 
level. Their masses are generated at one-loop by diagrams involving the sterile 
neutrinos (sj) on internal lines. The pattern of symmetry breaking induced by 
the scalar vacuum expectation values is 



SU(A) C x SU{2) L x SU{2) R x CP SU(A) C x SU{2) L x SU{2) 

^ SU(3) c xSU{2) L xU(l) Y 
^ SU(3) c xU(l) Q . 



(9) 



If we use as inputs at the electroweak scale a -1 = 128.5, = 0.23, and a^ 1 = 
8.33, and the boundary condition g^Vn) = ^(fj?), then using the one-loop (3 
functions with only gauge boson and fermion contributions, gives us Vr = 10 14 
GeV. The masses of the sterile neutrinos are required to be close to Vr in order 
to optimize the radiative mass generation. Thus we choose vq = 10 15 GeV. All of 
the scalars have a mass ~ Vr except for L\ and L\ which remain light and form 
the components of the Standard Model Higgs boson. 



3 Baryogenesis 

In this section we work in the symmetric phase of the Standard Model, i.e., 
vl = 0. We shall denote the light scalar doublet with components (X*, L\) as $*, 
the (complex conjugated) Standard Model Higgs boson. The out-of-equilibrium 
decay that we are then concerned with is that of the lightest (right-handed) 
SU{2) x U(l) singlet neutrino n\ into a left-handed lepton, l L , and the Higgs 
boson, $ |J. In order to calculate this, we first need to re-express the Yukawa 
interactions of Fig. 1(a) involving the singlet neutrinos Sj (i = 1,2,3) in terms 
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of the physical (Majorana) neutrino states rij (j = 1,...,6) obtained when one 
diagonalizes the neutrino mass matrix [Fig. 1(b)]. 




Fig. 1 Yukawa couplings of the heavy singlet neutrinos to a light left-handed lepton 
doublet and the Standard Model Higgs. The arrows indicate flow of left-handed isospin; 

and hij are coupling constants, (a) Coupling of the Sj. (b) Coupling of the phys- 
ical states rij. 



The vacuum expectation values of Eq. ([8]) along with the Yukawa couplings 
of Eq. ([7p lead to the following 6x6 mass matrix for the right-handed neutrinos 
in the basis (iVj, Sj) with i = 1, 2, 3: 





/ 

















2V2 




M = 


V 


(««) T uji 

2^2 


m ie iai 






mo2e lQ2 






m 03 e * a3 



where m 0i = yiVo/y/2. The matrix M will be diagonalized by a 6 x 6 unitary 
matrix U with M = UMe>U t where Mjj is diagonal, real and positive. 

We can now obtain the Yukawa couplings of the physical heavy neutrinos, n J ", 
with the light left-handed leptons, /*, and the Higgs boson $: 

C' Y = —h ij T\ j n j R $ + h.c 

where 

ha = E^ium+sj (12) 

k=l 

is now a complex 3x6 matrix, and rij are the 6 physical right-handed neutrinos 
arranged so that mi < m 2 < ... < m 6 . 

From here the calculation proceeds just as the cases presented in the literature 
(except for the fact that we have to account for 6 heavy neutrinos rather than 3), 



(11) 
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depending on whether we assume a period of inflation in the early universe [13 
or not H PA , and we now illustrate with a numerical example that the model 
can naturally generate a baryon asymmetry of the correct size. 

We choose as inputs the following matrices of Yukawa couplings (we will 
motivate this choice in the next section): 

/ 0.04 0.03 0.06 \ / -0.06 0.03 -0.04 \ 

4 = 0.06 0.42 0.24 ; k x r = 0.06 0.20 0.16 (13) 
V 0.06 0.08 3.5 / V 0.06 0.08 3.5 / 

y x = V2 x 0.025; y 2 = V2 x 0.05; y 3 = V2 x 0.35. (14) 
We make a simple choice for the phases «j of Eq. (||) : 

— 7r < a\ < 7r; a 2 = « 3 = 0. (15) 

Since rii is a Majorana fermion it can decay via the Yukawa interaction of 
Eq. ( pi) into both leptons and anti-leptons, i.e., 

n x -> Zi + $ 

- (IIT + ®*, (16) 
thus violating lepton number in its decays. The tree level width of ni is then 

r = (17) 

87T 

where h is defined in Eq. (|12D and 

(rfh) n = 3 x 10" 5 , (18) 

and 

mi = 4 x 10 11 Gel/ (19) 

for the inputs we have chosen. 

Further, if CP is violated there will be an asymmetry in the decay rates into 
leptons and anti-leptons. Let us define an asymmetry 

r _ t cp 



r + r CP 



(20) 



where T is the decay rate into leptons, and T CP into anti-leptons. Interference 
between the diagrams of Figs. 2(a) and 2(b) give rise to 8 which is calculated to 



be B 14 



$ = 2 ^ ( ^ )ii Elm[(^) lj ] 2 /K 2 Mi), (21) 
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where 



/(*) 



1 - (1 + x)ln 



1+x 



X 



(22) 



is proportional to the imaginary part of the loop amplitude of Fig. 2(b), and 
arises when the 4 and $ in the loop go on-shell. 

li 



(a) 



ni 




(b) 



n. 




Fig. 2. Decay of the lightest physical singlet neutrino n\ into a left-handed lepton and 
the Standard Model Higgs. (a) Tree level, (b) One-loop, with on-shell intermediate 
states. 



The calculation now splits up depending on whether we want to incorporate 
the effects of inflation or not. We will first estimate the baryon asymmetry in 
the inflationary case (although this calculation is more speculative than the non- 
inflationary one). We will then adapt the result to the case of no inflation, i 

Quantum fluctuations of the scalar field r\ (the inflaton), that drives inflation, 
give rise to density perturbations in the early universe which in turn lead to 
anisotropies in the cosmic microwave background radiation. Thus one can use 
the COBE measurements of the quadrupole moment of the microwave background 
T6| to constrain the parameters of inflation. For a generic inflationary potential 



as defined in Ref. ||13| , one finds 



Mt 



10" 



(23) 



where \x is the inflation scale, and Mp = 10 19 GeV is the Planck mass. Thus 
inflation takes place at a scale \i ~ 10 15 GeV. 

The fact that inflation occurs at or below Vq = 10 15 GeV, the scale at which 
CP is spontaneously broken, solves the problem of domain walls in this model. A 
domain wall arises whenever a discrete symmetry is spontaneously broken, and 



contains an unacceptably large energy density [17|. If one existed in our Hubble 



2 For a text book discussion of some of the estimates we make in this section and their 
accuracy, see Ref. |15|. 



7 



volume, it would have over-closed the universe many times over [I8|. Inflation at 
a scale below that of the discrete symmetry breakdown has the effect of diluting 
the density of domain walls to the level where it is extremely unlikely to find one 
in our observable universe. 



The mass of the inflaton is 



n 2 

mr] ~ f- ~ 10 n GeV (24) 



M ± 

and the reheat temperature is 

Trh * (r.Mp) 1 / 2 ~ ~ 10 7 - 10 8 GeV (25) 

where T v ~ m^/Mp is the inflaton decay rate, and T RH is defined as the equilib- 
rium temperature of the relativistic decay products of the inflaton. 

Thus m v ~ 10 11 GeV >, mi, and the inflaton is heavy enough to decay into the 
right-handed neutrino n\ whose subsequent decay generates a lepton asymmetry. 
A simple way to estimate the baryon asymmetry in this case is |13| 







Pn s 


s 


s 


m v s 




Trh , 

~ 


~ (10" 









10~ 4 )S. (26) 

Here we have used the approximation that all of inflaton's energy density, p V) is 
instantaneously converted into the energy density of relativistic particles: p v ~ 
g*T RH . The subsequent entropy density is s ~ 9*T R [ H (g* is the number of 
relativistic degrees of freedom), and m v and T RH are obtained from Eq. (^4j) and 
Eq. (p5|) above. The baryon asymmetry is then 

% 3 = -^x7^~ -(10- 3 - 10- 4 )5 (27) 

where we have used 



n 7 79 s 



ub_ = 28 n [B -L) , , 

s 79 s K ' 

as a result of the electroweak B + L violation |19| , and the fact that s = 7ra 7 



today. In Fig. 3 we plot the CP asymmetry 5 as a function of the phase oc\ for 
the choice of inputs of Eqs. (fl3| , |i~4j) . We see that the choice ai = — it/ 2 gives 



us 5 = —3 x lO -7 and subsequently tjb — 10 -10 in agreement with the observed 
value. We pick the value with |ai| = 7r/2 as it might be consistent with the 
notion of a "maximal" CP violation in the early universe. 
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Fig. 3. The CP asymmetry 8 vs. the phase a L 
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The last thing we need to ensure is that the lepton asymmetry is not washed 
out by the lepton number violating process IlIl — > of Fig. 4. 




n. 




hS X 1l 

Fig. 4. Induced AL = 2 interaction. The arrows indicate flow of lepton number. 



This requires that the process be out of equilibrium when the electroweak 
B + L violation is in equilibrium ]7|, [13], |I§fl , which happens at a temperature 
Tr ~ a 2 T^H ~ 10 5 GeV, where a is a typical gauge coupling |p0| . The bound 
one needs to satisfy is 



> 1.4 x 10^ \/T R Mp (29) 



ii 



which is easily satisfied for the parameters we use. 

Although we prefer the inflationary scenario because it cures the domain wall 
problem as mentioned above, and because of its other attractive features, we now 
show how to adapt the calculation to the case of no inflation. 

In the non-inflationary case, we first need to check if n\ is out of equilibrium 
when it decays i.e. 

T t < H(T = mi) (30) 



where H is the Hubble constant. This translates to 

(tfh) u m 1 < 20ml 
8^ ~ M P 



(31) 



which is satisfied for our choice of inputs. The lepton asymmetry is then given 
by 



til n^, 

n*J L A r^j <__ 

-2 



10~ 2 5 (32) 
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where we have used g* ~ 100 between the scales vl and vr. Thus we see that 
in this case a sufficient baryon asymmetry can be achieved with a smaller CP 
violation than in the inflationary case. This is borne out by Fig. 3 where we 
see that smaller values of a.\ indeed lead to smaller values of 5. If we choose 
aii — — tt/96, we get 5 = — 3 x 10 -9 and r] B ~ 10~ 10 . 

Once again we need to ensure against a wash-out of this symmetry. In this 
case, the B + L violation is in equilibrium up to a temperature of about 10 12 
GeV, and hence in the RHS of Eq. ( p9|) we use the temperature at which the 
lepton number violation occurs: Tl = m\ = 10 11 GeV. The inequality 



m l ^ i „ tn-2 



ii 



> 1.4 x 10~ 2 JT L M P (33) 



is indeed satisfied, and once again we find that our model can generate a sufficient 
baryon asymmetry. 



4 Fermion Masses and the CKM matrix 

The calculation for the fermion masses and mixings exactly follows that in |§, 
except that there are now 3 generations of scalars rather than 2, and that we now 
allow for CP violation, hence the orthogonal matrix O in the mass formulas is 
replaced by the unitary matrix U, or U* as appropriate. We relegate the details 
to Appendix B, restricting the discussion in this section to qualitative features, 
and demonstrating with a numerical example. 

As mentioned earlier, none of the Standard Model fermions get masses at 
tree level, even after the standard electroweak group SU{2) L x U(l)y is broken. 
However, one-loop masses are generated by the diagrams of Figs. (5). X in 
Fig. 5(a) is an SU(4) lepto-quark gauge boson, and Z' of Figs. 5(e), 5(f) is 
the neutral gauge boson that couples to the linear combination of the broken 
diagonal generators U{V)b-l an d U(1)r of SU(4) and SU{2)r respectively, that 
is orthogonal to the hypercharge U{1)y- Z of Fig. 5(f) is just the Standard 
Model Z boson. 

One can estimate the magnitude of the mass terms from Figs. 5(a) and 5(e) 
involving gauge boson exchange as 

9\92 mo 9x92 , a\ 
(An) 2 rriQ (Airy 

where g\ and are the gauge couplings at the two vertices. The factor of (47r) 2 
in the denominator is from the loop integral, and we have replaced the sterile 
neutrino mass tuq by Vr. 
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Fig. 5. One-loop processes that generate the fermion masses. The arrows here indi- 
cate the flow of left or right-handed isospin. (a) Up-type quarks with gauge bosons, (b) 
Up-type quarks with scalars. (c) Down-type quarks, (d) Charged leptons. (e) Dirac 
mass for neutrinos, (f). Majorana mass for neutrinos. 
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Notice that if we replace vl on the external leg of Fig. 5(e) by the unshifted 
Higgs field, $, we get an effective interaction similar to that of Fig. 1(a), but with 
Nj on the external leg. One might then worry about whether it was correct to 
ignore this contribution to the baryogenesis calculation of the previous section. To 
estimate the magnitude of the effective Yukawa couplings of the Nj to the k and 
$, one can simply divide the result of Eq. ( |34"D by v^. If we use gig 2 /(4:7i) ~ 1/40, 
we get effective Yukawa couplings y ~ 0.002 which are smaller than the couplings 
k\ of Eq. (|T3|) by an order of magnitude. Thus their contribution to the total 
width of ri\ can be be ignored. The fact that this process makes no contribution 
to the CP asymmetry can be understood by replacing all the particles in Fig. 5(e) 
by physical particles. Then, with an on-shell n\ on the external leg, and rii and 
Z' in the loop, the diagram can never have an imaginary part since ri\ is the 
lightest physical neutrino, and the particles in the loop cannot be on-shell. These 
estimates were confirmed by a direct calculation of the appropriate three-point 
form factors, and their contributions to the baryogenesis calculation. 

Using arguments similar to the ones for Figs. 5(a) and 5(e), one can estimate 
Fig. 5(f) to give 

m G2 ~ y~j Y2 — v L , (35) 
(Air) 2 Vr 

and the scalar exchange contribution of Figs. 5(b), 5(c) and 5(d) to the fermion 
masses are approximately given by 

3 a K a 

One immediately notices a fundamental difference between fermion masses 
in this model and those in the Standard Model. Here all the fermion masses 
are proportional to the squares of coupling constants as opposed to depending 
linearly on them.! This has the effect of reducing the hierarchy in coupling 
constants needed to reproduce the large hierarchy in the observed masses. 

If this model were supersymmetric then the estimates for the fermion masses of 
the previous equations would be scaled by a factor of Msusy I vr for Msusy < v r, 
where Msusy is the supersymmetry breaking scale. Thus the discovery of TeV 
scale supersymmetry would rule out this model for generating the fermion masses. 

Although all the quarks and leptons of one generation have identical Yukawa 
couplings in this model, the scalars to which they couple have different masses 



3 A more careful estimate would show that the gauge contribution of Eq. (|3J) is numerically 
similar to the scalar contribution with k\ and k} r of Eq. (|3^) . This must be since some of the 
scalars L 1 and R 1 form the longitudinal components of the massive gauge bosons, cf. Eqs. (^,[|). 
The choice of "gauge" or "scalar" contribution then depends on the gauge we calculate in, with 
the final answer being gauge independent. 
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and mixings determined by the parameters of the scalar potential (the explicit 
formulas for the masses are listed in Appendix B). Using this fact, and some of 
the features of Figs. 5, one can qualitatively understand how this model gives rise 
to the rather complicated observed spectrum of fermion masses. 

The up-type quarks get masses from both gauge boson and scalar exchange 
[Figs. 5(a), 5(b)]. This larger number of diagrams, and the possibility of con- 
structive or destructive interference between them allows us to generate the large 
hierarchy in their masses. In particular all the diagrams that contribute to the top 
quark mass interfere constructively to generate the large observed value, whereas 
destructive interference between various diagrams is responsible for making the 
up quark lighter than the down quark. 

The down-type quarks and charged leptons both get masses only from scalar 
exchange, explaining the similarity in their spectra [Figs. 5(c), 5(d)]. 

The neutrinos get masses from gauge boson exchange [Figs. 5(e), 5(f)]. How- 
ever, the large mass of the sterile neutrinos results in see-saw suppression of the 
physical left-handed neutrino mass to a level compatible with the MSW plj so- 
lution to the solar neutrino problem. The neutrinos could in principle also get 
masses from scalar exchange, but this contribution vanishes for a natural choice 
of parameters in the scalar potential. 

The complex phase in the CKM matrix arises as a result of the single phase 
«i in the sterile neutrino mass matrix. 

Let us illustrate these features with a numerical example. We choose the 
following Yukawa couplings in addition to k\ 



K 



R- 



Vi of Eqs. (|13|, H 





/ 0.04 


0.03 


0.06 


4 = 


0.03 


0.4 


0.2 




V o.i 


3.5 


3.5 




1 -0.2 


0.02 


-0.2 


4 = 


-0.2 


0.4 


-0.4 




v 3.5 


3.5 


3.5 



K 



R 



0.04 


0.03 


0.06 


-0.03 


0.4 


-0.2 


0.1 


3.5 


3.5 


0.2 


-0.02 


0.2 


-0.4 


-2.0 


-1.0 


3.5 


3.5 


3.5 



(37) 



The procedure we used to pick the Yukawa couplings was the following. The 
diagonal entries of # were arbitrarily chosen to reflect the hierarchy 1 : A : A 2 
with A = 0.1, going from the third generation to the first, in accordance with our 
policy of not allowing too large a hierarchy in the couplings. The scale was set by 
the (3, 3) entry which was picked to be 3.5, the largest value with which we felt 
comfortable doing perturbation theory. The first row and column were chosen to 
reflect the magnitude of the (1,1) element, and the (2,3) and (3,2) elements were 
arbitrarily chosen to reflect either the (2,2) or the (3,3) elements. Some of the 
entries were randomly assigned minus signs (except the entries in the third row, 
which contribute to the top quark mass, and were all required to be positive). For 
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the k\ r , all the elements of the third row were set to 3.5 in order to maximize 
the top quark mass, and the elements of the first and second rows were varied in 
order to fit the data. 

Given the inputs above, for a choice of parameters of the scalar potential 
similar to those of one obtains the following fermion masses and mixings (at 
the electroweak scale ~ 100 GeV). 

m u = 7.0 MeV; m c = 1.2 GeV; m t = 160 GeV (38) 

m d = 8.0 MeV; m s = 140 MeV; m b = 3.9 GeV (39) 
The absolute values of the CKM matrix elements are 

0.98 0.21 0.004 
\V\ = | 0.21 0.98 0.06 j (40) 
0.01 0.06 1.0 

As a measure of the CP violation in the CKM matrix we first calculate the 
parameter Jcp, where 

Jcp = lm[V cb V us V: b V*} (41) 

is twice the area of the unitarity triangle of CKM matrix elements, and is a rephas- 
ing invariant measure of the CP violation in the CKM matrix with 3 generations 
of quarks. We plot Jqp as a function of the phase ol\ in Fig. 6. In the Wolfenstein 



parametrization of the CKM matrix [p2 |, we have Jcp = A 2 X 6 r] ~ 4 x 10 5 where 



we have used A = 0.85, A = 0.22, and rj ~ 0.5. For our choice of inputs, and 
with cti = — 7r/2 we see from Fig. 6 that 

J CP = 3 x 10" 5 , (42) 

which has the correct magnitude and sign. 

Notice that if we attach a photon to the charged particles on the internal 
lines of the diagrams in Fig. 5, we could induce an electric dipole moment for 
the quarks (and charged leptons) at one-loop order. This is to be contrasted 
with the Standard Model, where the fermion electric dipole moments arise only 
at three- loop order p3, 24]. An estimate of the dipole moment is then 



dii> ~ -~t~ e ' cm (43) 

where is the mass of the fermion in question. This is well below the estimated 
Standard Model value for the quarks fl25 |, and hence the low energy hadronic CP 



violation in this model is effectively of the CKM type. 
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In order to estimate the kaon CP violating parameter e, we assume that the 
top quark contribution dominates the box diagram that leads to K° — K° mixing. 
In this case one has fEEL E7H 



C e B K J C p\V ts \ 
~3x 1(T 3 



where we have used 



a 



G 2 F F 2 K m K Mlr 
6V2n 2 AM 



3.8 x 10 4 , 



(44) 



(45) 



and Bk = 0.65 [Bk measures the accuracy of the vacuum saturation approxima- 
tion) . 

In the lepton sector we obtain 



0.7 MeV: 



m, 



90 MeV: m T 



7 x 1(T 5 eV; m u = 4 x 10" 3 eV; 



m,. 



The absolute values of the lepton mixing matrix are 



Wv 



I 0.92 0.39 0.03 
0.39 0.91 0.12 
V 0.03 0.12 0.99 



1.3 GeV. (46) 
= 3 x 10" 2 eV (47) 

(48) 



We should point out that since the neutrino masses depend sensitively on the scale 
vr via the see-saw mechanism, the above values should be taken as representative 
of the range the actual masses should lie in. This prediction of the scale of the 
neutrino masses is then fairly robust since they are generated only by gauge boson 
exchange, and hence independent of the Yukawa matrices Kl'r and details of the 
scalar potential (they do depend on k} l r which, however, are constrained by the 
baryogenesis calculation of Sec. 3). It is interesting to note, then, that the v e 
and masses lie in the correct range to explain the solar neutrino deficit as a 
result of the MSW effect, with the mixing angle being compatible with the 



large angle solution p8 



We should point out, however, that the prediction for 
this mixing angle is not as robust as that of the masses since it is dominated 
by the contribution of the unitary matrix that diagonalizes the charged lepton 
mass matrix, and hence dependent on details of the scalar sector. Although the 
Vjx — v T mass difference lies in the correct range for the proposed solution to 
the atmospheric neutrino problem, the required mixing angle is generically much 
smaller than that required by the theoretical fits to the data J28| . 
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This model also predicts CP violation in the lepton sector, due to complex 
phases in the lepton mixing matrix. However, these effects are suppressed by the 
smallness of the neutrino masses, and hence, unobservable. If one applies the 
estimate of Eq. ( [43|) for the electric dipole momments of the charged leptons, one 
gets a value larger than what one would calculate just using the phases in the 



mixing matrix ||23|| , however it is still much to small to be observed. 

Thus we are able to generate realistic fermion masses and mixings within our 
stated objective of working in a quark-lepton symmetric model without a large 
hierarchy of coupling constants. Not only that, but starting with the same phase 
ai = — tt/2 in the sterile neutrino mass matrix that generates a lepton (and 
hence baryon) asymmetry, we are also able to generate the correct amount of CP 
violation in the CKM matrix Eq. (pE2"|). 

All of the inputs to the model were defined at the high energy scale vr = 10 14 
GeV. We would like to remark on the procedure we have followed in obtaining 
the outputs at 100 GeV. The effective theory below the scale vr is just the 
Standard Model and hence we run the masses using Standard Model (3 functions. 
Effectively this results in a scaling of the light quark masses by a factor of 2.5 
to account for the QCD effects, whereas the top quark mass scales by a factor 
of 1.7, on account of its slower running due to its large Yukawa coupling to the 
Standard Model Higgs. None of the mixing angles, nor the lepton masses run to 
the degree of accuracy we're concerned with. 

As we have mentioned, the discovery of low energy supersymmetry would 
rule out this mechanism of mass generation. However the scenario whereby a CP 
violating mass matrix for sterile neutrinos is the source of both the baryon asym- 
metry and the CKM phase may still remain viable. We are currently studying 
the possibility of building this scenario into a supersymmetric model. Another 
possibility we are studying is the extension of this model to include family sym- 
metries for the Standard Model fermions in order to reduce the arbitrariness of 
the Yukawa couplings. The fact that the Yukawa couplings in this model only 
range from 0(1) to 0(X 2 ) with A = 0.1, may allow a simple assignment of fam- 
ily charge to the Standard Model fermions, and the possible use of non-abelian 
family symmetries P5fl. 



5 Conclusions 

We have presented a simple extension of the Standard Model based on the Pati- 
Salam gauge group SU (4) x SU{2) l x SU (2)r. CP is spontaneously broken in the 
early universe, and the CP violation manifests itself as a single phase of —ir/2 in 
the mass matrix of the sterile neutrinos that are present in our model. With this 
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phase as the source, we can generate both the baryon asymmetry of the universe, 
t]b — 10 -10 , and e ~ 10~ 3 , the CP violating asymmetry in kaon decays. 

Masses for the fermions are generated radiatively at one-loop order, and we 
can obtain realistic masses and mixings without the large hierarchy of coupling 
constants required in the Standard Model. The neutrino masses lie in the correct 
range for the MSW solution to the solar neutrino problem. There will be CKM- 
type CP violation in the lepton sector, however the smallness of the neutrino 
masses makes any effects unobservable. 
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Appendix A: The Scalar Potential 

Given the hierarchy t> > Vr ^> vl, the simplest way to study the scalar potential 
is to analyze the three sectors comprising a 1 , R l , and L l independently. While 
this suffers from a lack of generality, it makes the analysis extremely transparent 
without affecting our results. 

For the a % we work with a simplified version of the Z3 invariant scalar potential 



ofRef. 12 



V(a) = -2m?(<7l<7i) - 1m\{a\a 2 ) - 2m\{(j\^) 

+a 1 (a\a l ) 2 + a 2 [a\o 2 f + a z [o\a z f 
+ [C 1 (a 1 V 2 )( ( T 1 V 3 ) + C 2 (o\o x )(o\oz) + Ci(a\a^(a\a % ) + h.c] (49) 

We then assume the vacuum expectation value (a)j = Vje tai / 'a/2, and study the 
equations to minimize the potential. In particular, for the choice C 2 = C3 = 
0, rrii ~ 10 15 GeV, we recover 

vl = v 2 2 =vl = vl = x — -; oty = --, a 2 = a 3 = (50) 

which are the values used in Eqs. |i~5D . 
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For the R l and L l , we specialize to the one-family case, since L 1 = L and 
R 1 = R are the only scalars that get vacuum expectation values (Eq. (§)). 

V(i?) = —2ii 2 R R ia R ia + X R1 R ia R ia RjpR^ + X R2 R ia Ri a R^ Rjp + ^m^aR^RpRj 

(51) 

= -2fi 2 L L ia V a + X L1 L ia L ia L jf3 L^ + X L2 L ia L ja L ifi L jl3 + X L3 L ia L ja V p L? . 

(52) 

Here i = 1,2 is the SU(2) L or 5*7(2)^ index, and a = 1,2,3,4 is the 517(4) 
index. L lQ = (L ia )* and = e l ^Lj a , and similarly for L — > R. We can then 
generalize the arguments of Ref. [S(| to show that if the following conditions are 



satisfied: 

A,R2 < 0; \m + Ar 2 > 0; Xr 3 > Xr 2 or \^m\ > 2A^x + Xr 2 
X L 2 <0; Ali + A L2 >0; A L3 > A L2 or \X L3 \ > 2X L1 + X L2 (53) 
the absolute minimum of the potential is at 

/ p\ ( v R /V2 \ 

{R) = [o ) 

(L) = ( o o o o J (54) 

with 

(Am + X R2 )v 2 R = 2/4 

(A L i + A L2 )t;| = 2/il. (55) 

For fj, R ~ 10 14 GeV and Hl — 100 GeV, we then obtain v R ~ 10 14 GeV and 
vl — 100 GeV. There will be other constraints on the parameters of the scalar 
potential to ensure (L)- 7 = (R)i = for j ^ 1, which are easily satisfied for 
natural choices of the parameters. The only light scalars are then L\ and L 4 , 
with the rest having masses ~ v R . 

We would like to point out that in order for scalar exchange to generate 
fermion masses one actually needs L % — R l mixing to interchange left and right 
weak isospin as one can see from Figs. 5(b), 5(c), 5(d). What we do in the actual 
calculation is to allow individual couplings to be non-zero, but require that the 
combination of constants that couple vl to v R in the equations that determine 
the minimum to add up to zero, thus leaving Eq. fl55|) unchanged even in the 
more general case (see Ref. || for details). 

Another point we would like to make is that one could in principle start with 
a scalar potential that was L — R symmetric, and use the couplings of the L l , R 1 
to the a % to break this symmetry as in Ref. [BT[ . 
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Appendix B: Details on the Fermion Masses 



In this appendix we would like to fill in some of the details, and list the explicit 
formulas used to calculate the fermion masses. 

All the calculations were done at zero external momentum, and with physical 
particles propagating in the loop diagrams of Fig. (5). Thus, before we get to the 
formulas for calculating the Standard Model fermion masses, we need to obtain 
the masses and mixing angles for the singlet neutrinos, gauge bosons, and scalars 
that propagate in the the diagrams of Fig. (5). 

The vacuum expectation values for the scalars, a 1 , L 1 , and R 1 [Eqs. 



5C 



along with the Yukawa couplings of Eq. (0) generate the following 9x9 

for % 



mass matrix for the neutrinos in the basis (z/ c , iVj, s 
( 



1,2,3, 



M 



\ 









(4H 
2V2 








2V2 


2V2 


2V2 


m ie iai 






mo2e"* 2 






mo3e iC * 3 



(56) 



where m 0i = yiV /y/2. The matrix M will be diagonalized by a 9 x 9 unitary 
matrix U with M = UMe>U t where is diagonal, real and positive. This 
matrix has eigenvalues mi = rri2 = = corresponding to the physical states 
ni, n>2, n-3 that are mostly admixtures of the Vi, and Tnt±,...,m§ ~ vr for the 
physical states 714, Tig that are mostly admixtures of the Ni and Sj. The mixing 
between the (ui) and the (Ni, Sj) is of order vl/vr. 

The gauge bosons for the gauge group SU(A) x SU(2)l x SU{2)r are: 



G„ 



G 



8/1 



V2G 



+ 



V6 



V2G 



12/i 



V2G 



13/i 



12/1 

V2G 



+ 



V6 



V2G 

V2G 
2G, 



13/t 



V2X 



i/i 



\ 



23/t 



23/1 



+ 



B„ 



V2X 



l/i 



y/2X. 



2/1 



V2X- 



3/1 



(57) 



V6 / 



W, 



V2W, 



V2W, 



and 



V2W, 



fi/i 



Lfj, 



(58) 



(59) 



23 



where the are the gluons, 



is the diagonal gauge boson that couples to 



L, and the are the lepto-quarks. 



B 

As a result of the spontaneous symmetry breaking outlined in Eq. 
charged gauge bosons get the following masses: 



M* x = 9 -f[v 2 R + vl] 



M 



^ L~tt ' ~U\1 —Wl ^ ~L,1 —WR | 

The neutral gauge bosons have the mass squared matrix 



£L,2. M 2 

■ V Li Wr 



9R2 

J li- 



the 



(60) 



/ 



o 



o 

9rVr 



-(3/V6)g L g s v 2 L 
-(3/V6)g R g s v R 



(61) 



V -{3/VB)g L gsvl -(3/Ve)g R9s vl (3/2)g 



in the basis (W®, W R , B°), with eigenvalues 



M 2 = 0; Mf 



W ? r 9 



2^ + 3^1 



■O(^K)]; Mf, = f[2g R +3gl + 0(vl/v R )} 

(62) 



and eigenvectors 



7 




—twSw —tw\ / c 2 w 
—\/czw/ cw tw 



wl 

W° R 

Br, 



(63) 



For vl <C v r , the usual electroweak relation M w 
where we define 



M 2 r 2 
IV1 Z C W 



2 _ c 2 



3^1 



3^1 + igisrk + Wig's 



2 „2 



,2 „2 



is still maintained 



(64) 



For the scalars we make the simplifying assumption that the different gener- 
ations don't mix. In this case, within each generation, we will get a 2 x 2 mass 
matrix for each of the differently charged scalars. As an example of the notation 
we will use, consider the scalars L\, R\. Diagonalizing their mass matrix will give 
us two physical states with masses M u3 , M' u2i and a mixing angle s u3 . The masses 
will generally be of order v R and the mixing angles of order Vl/v r . Care should 
be taken to spot the scalars that get absorbed to form the longitudinal compo- 
nents of the massive gauge bosons, and to treat them appropriately. Since we do 
our calculation in 't Hooft-Feynman gauge, these scalars are simply assigned the 
mass of the corresponding gauge bosons. 

Having enumerated the physical neutrinos, gauge bosons, and scalars, we can 
finally list the formulas used to calculate the fermion masses (the indices (a, b) 
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in the following formulas are generation indices). For the up-type quarks we get 
from Fig. 5(a): 

9 9 9 9 

(7 _ 771 771 

"•"'"^S^'Wn.I^rM^)]. (65) 

From Fig. 5(b) we get: 



rn 



3 ^ara ^.bn 
ab \- ^LjJ^Rj 



m,n=l V^" / j=l 

M 2 . m 2 M' 2 m 2 

"J lnf— — — ) — ^ ln(— Ml (66) 



^f-M 2 , V M 2 / m 2 -M^. V M^ 



and we sum over j — 1, 2, 3 for the three generations of scalars. 
For the down- type quarks, Fig. 5(c) gives us 



3 ^am „bn 9 



m,n=l V 47r J i=l 



Mj- m 2 m 2 

[ IZ* ln(^) — ^ ln(^)l (67) 

[ m 2 -M 2 , mi M 2 / m?-M'r n[ M'J )l - 



dj i dj dj 

For the charged leptons, from Fig. 5(d), we get 

3 „am ^.bn 9 



ab v ^—3Lc .vrr r/* ™ 

'./ — b ej 2-^ U m+6,i u n+6,i' ri t 



m,n=l V*" ; i=l 



M 2 . m 2 M' 2 m 2 

Here, however, there is no contribution for j = 1 since the L\ and R\ form the 
longitudinal components of the Wl and Wr, respectively, and don't mix at tree 
level. 

For the neutrinos, Fig. 5(e) gives us the following Dirac mass matrix 
where 



3gs 3t w gs . Ql\Jciw ,„ n \ 

Jl = 76 ; 92 = ^T + ^— (70) 
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Fig. 5(f) leads to the following Majorana mass matrix for the left-handed neutri- 



nos, 



m ab =2 7 ^E^ i m i [- f ^ln(^)] 
i=i m i lvl z m : 



(4tt 



where 



+2 (Js^A M 3 )] (71) 

9o = — • (72) 
c w 

There will be a similar Majorana mass matrix for the right-handed neutrinos, 
N a , with gi — > g 2 , and a, b — > a + 3, b + 3 in the first term of Eq. (|7lD , and no 
contribution from the second term. These 3x3 mass matrices for the neutrinos 
fill in the appropriate zero blocks in the 9x9 mass matrix of Eq. (|56|), which is 
then rediagonalized, resulting in small but finite see-saw suppressed masses for 
the left-handed neutrinos. 
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